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We calculate high-order quasinormal modes with large imaginary frequencies for electromagnetic 
and gravitational perturbations in nearly extremal Schwarzschild-de Sitter spacetimes. Our results 
show that for low-order quasinormal modes, the analytical approximation formula in the extremal 
limit derived by Cardoso and Lemos is a quite good approximation for the quasinormal frequencies 
as long as the model parameter riKi is small enough, where ri and ki are the black hole horizon ra- 
dius and the surface gravity, respectively. For high-order quasinormal modes, to which corresponds 
quasinormal frequencies with large imaginary parts, on the other hand, this formula becomes inac- 
curate even for small values of riKi. We also find that the real parts of the quasinormal frequencies 
have oscillating behaviors in the limit of highly damped modes, which are similar to those observed 
in the case of a Reissner-Nordstrom black hole. The amplitude of oscillating Ke{uj) as a function 
of Im(cij) approaches a non-zero constant value for gravitational perturbations and zero for elec- 
tromagnetic perturbations in the limit of highly damped modes, where oj denotes the quasinormal 
frequency. This means that for gravitational perturbations, the real part of quasinormal modes of 
the nearly extremal Schwarzschild-de Sitter spacetime appears not to approach any constant value 
in the limit of highly damped modes. On the other hand, for electromagnetic perturbations, the 
real part of frequency seems to go to zero in the limit. 
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I. INTRODUCTION 

The quasinormal modes (QNMs) of spacetimes con- 
taining black holes have been studied since the pioneering 
work of Vishveshwara Q , who first observed quasinormal 
ringing of a Schwarzschild spacetime in his numerical cal- 
culations. The main motivation to study the QNM of a 
black hole is twofold: One is to answer the question of 
whether the spacetime is stable or not, and the other 
to know what kind of oscillations will be excited in the 
spacetime when some perturbations are given. In an as- 
trophysical point of view, the latter is quite important 
from the observational point of view because we could 
determine fundamental parameters of a black hole, such 
as the mass or the angular momentum, through the in- 
formation of the QNMs. Thus, a large number of studies 
on the QNMs of spacetimes containing black holes have 
been done (for review, see, e.g., P,^). 

There is another interesting aspect of the QNMs in 
a black hole which is related to the quantum theory of 
gravity. Bekenstcin and Mukhanov discussed the rela- 
tionship between the fundamental area unit in the quan- 
tum theory of gravity and a Bohr transition frequency, 
applying Bohr's correspondence principle with a hydro- 



gen atom to the quantum theory of a black hole 0, |^ . 
For a Schwarzschild black hole, they then derived the 
Bohr transition frequency u, given by 



LO = Ink/SirM , 



(1) 
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where k — 2,3,4- •■, and M stands for the mass of a 
black hole. Furthermore they predicted the value of k to 
he k — 2, and suggested that this frequency should be 
equal to classical oscillation frequencies of the black hole. 
A few years ago. Hod however noticed that ii k — 3, the 
frequency given by formula JQ) is in quite good agree- 
ment with the asymptotic frequency of the QNM of a 
Schwarzschild black hole in the limit of highly damped 
modes, and proposed to apply Bohr's correspondence 
principle in order to determine the value of the fundamen- 
tal area unit in the quantum theory of gravity, namely 
the value of A: (see, also 0). Since Hod's proposal, 
the QNMs with large imaginary frequencies of spacetimes 
including black holes have been attracted much atten- 
tion, and a lot of papers related to this subject have 
appeared in order to see whether Hod's conjecture is ap- 
plicable not only for a Schwarzschild black hole but also 
for other black hole spacetimes. For example, Motl an- 
alytically obtained an asymptotic constant value of the 
QNM frequencies of a Schwarzschild black hole ^ , which 
had been obtained numer ically by Nollert Motl and 
Neitzke fl^, Berti and Kokkotas 0|, Neitzke 0|, and 
Andersson and Howls 13] studied the asymptotic behav- 
iors of the QNMs in a Reissner-Nordstrom black hole. 
Berti et al. Q, and Hod jl^ discussed the QNMs of a 
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Kerr black hole in the limit of highly damped modes. As 
for a Schwarzschild-de Sitter (SdS) black hole, Cardoso 
and Lemos jl^ and Maassen van den Brink analyti- 
cally obtained asymptotic form of the QNM frequencies 
in almost the extremal limit, in which the cosmological 
horizon becomes very close to the black hole horizon (for 
the case of a Schwarzschild black hole in an anti-de Sitter 
spacetime, see, e.g., 18, 19j). 



A large number of papers related to Hod's conjecture, 
which have recently appeared, suggest that Hod's con- 
jecture is not universal, at least as it stands, even though 
it is applicable for Schwarzschild black holes in four and 
higher dimensions (for the higher dimensional case, see, 
e.g., 10, 12, 27, 28J). For example, a real part of the 
quasinormal frequencies in a Reissner-Nordstrom black 
hole appears not to go to any constant value in the limit 
of highly damped mode, but shows some periodic behav- 
iors as the imaginary part is increased B IH El [11 . 
This means that Hod's conjecture is not applicable in the 
Reissner-Nordstrom case because a real part of the quasi- 
normal frequencies does not have a limit as the imaginary 
part goes to infinity. In such a situation, it is necessary to 
explore a problem whether or not there is another black 
hole spacetime in which Hod's conjecture is applicable. 
The purpose of this paper is to improve our understand- 
ing of this problem and we are concerned here with the 
QNMs of SdS spacetimes. SdS spacetime has no spatial 
infinity but has cosmological horizon, and if Hod's con- 
jecture is applicable, it is interesting to see whether Hod's 
conjecture depends only on the black hole horizon, but 
not the cosmological horizon. Furthermore, recent obser- 
vations show that the universe does have a non-zero pos- 
itive cosmological constant. Therefore, SdS spacetimes 
are considered to be a good simple model of a black hole 
in the universe. 



In this study, in particular, we calculate numerically 
the QNMs of nearly extremal SdS spacetimes for rea- 
son we describe below. We therefore assume the surface 
gravity ni at the black hole horizon to be ki < 10^^ r^^ , 
where ri stands for the coordinate radius of the black 
hole horizon. Note that the extremal limit of the SdS 
space corresponds to the limit of ki ^ 0. In a nearly 
extremal SdS black hole, as mentioned before, an analyt- 
ical formula for quasinormal frequencies can be derived 
[TgL IT^ . One of the aims of this paper is to examine 
whether this analytical formula is correct in the limit of 
highly damped modes. 



The paper is organized as follows. In §2 we present the 
basic equations for obtaining QNMs in the SdS spacetime 
using Leaver's continued fraction technique [20j, which 
was extended to the case of the SdS spacetime by Moss 
and Norman |2lj| . Numerical results are given in §3, and 
S4 is devoted for conclusion. 



II. METHOD OF SOLUTIONS 

In order to examine the QNMs of the SdS spacetime, 
we make use of the same formalism as that derived by 
Moss and Norman [23|, who obtained low-order quasi- 
normal frequencies of the SdS spacetime for a wide range 
of the model parameter. The line element of the SdS 
spacetime is given by 



ds'^ = -^dt'^ + ^dr'^ 



Hd0^ 



sm 



where 



A = - 2Mr - -Ar^ . 



(2) 



(3) 



Here, M and A stand for the mass of the black hole 
and the cosmological constant, respectively. If a non- 
negative cosmological constant is assumed, namely de 
Sitter spacetime, there are two horizons, whose radial 
coordinates are given as positive solutions of r^^A = 
r — 2M — ^Ar^ = 0, one is the black hole horizon r = ri 
and the other the cosmological horizon r = r2, where 
r2 > ri > 0. Note that the equation r - 2M - ^Ar^ = 
has one negative solution r — < for the SdS space. 
One of the important model parameters of the SdS space- 
time is the surface gravity ki, defined by 



1 dA 
Ki = hm TT^-r- ■ 
r^ri Zr"^ dr 



(4) 



In terms of the non-dimensional surface gravity riKi, the 
mass and the cosmological constant can be written as 



M 



1 



ninm + l), A = rf2(i_2riAti), (5) 



which shows that < riKi < 1/2 for the SdS spacetime. 
In this study, we employ the non-dimensional parameter 
riKi to specify the SdS spacetime. 

By virtue of the symmetry properties of the SdS space- 
time, the master equations for the scalar (s = 0), electro- 
magnetic (s = 1) and gravitational perturbations (s = 2) 
can be cast into a wave equation of the simple form, given 

by miilil 



d20(r) 



+ [uj^-V{rmr)=Q, 



(6) 



where denotes the tortoise coordinate, defined by 
dr^ = r^A^^dr, and lj is the oscillation frequency of the 
perturbations. Depending on the type of perturbations, 
here, the effective potential is explicitly given by 



2M 



1(1 + 1) ,. 



2Ar 
3 



for s = 
for 5 = 1 
for s = 2 , 



(7) 



where I means the angular quantum number of pertur- 
bations. Here, only the axial parity perturbations have 
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been considered for the gravitational case because quasi- 
normal frequencies of the polar parity perturbations are 
the same as those of the axial parity perturbations (for 
the proof, see Appendix). 

The QNMs of the SdS spacetime are characterized by 
the boundary conditions of incoming waves at the black 
hole horizon and outgoing waves at the cosmological hori- 
zon, given by 

0(r)-|Cr (8) 
' [ e'"'^' as r, — > — oo , ^ ' 

where the time dependence of perturbations has been as- 
sumed to be e*"*. In general, it is impossible to adapt 
the boundary condition ((HI in a straightforward numeri- 
cal integration to obtain quasinornial frequencies. Some 
special technique is therefore required for computations 
of QNMs. In the present investigation, we employ a stan- 
dard technique devised by Leaver, namely the continued 
fraction method [Mill- 

To apply the continued fraction method to the SdS 
spacetime, it is convenient to introduce a new indepen- 
dent variable, defined by x = . With this new variable 
X, the asymptotic form of the perturbations as r* — > oo 
can be rewritten as 

e"*"'^* = {x - x^yp^x - X2yHx - ^3)'"" , (9) 
where Xi = r^^ and pi ~ iu/ {2Ki) for i = 1, 2, 3, where 

Ki = M(xi - X2)(a;i - X3) , 

H2 = M(x2 - xi)(a;2 - xs) , (10) 

K3 = M(x3 - xi)(a;3 - X2) . 



The perturbation function is expanded around the 
black hole horizon as 




where — 1 and a„'s for n > 1 are determined by the 
three term recurrence relation, given by 



Q;„a„+i -I- /3„a„ -I- 7„a„_i , (12) 

where 

an = 2M(xi-xz){n'^ + 2{pi + l)n + 2pi + l] ili) 

fin = -2Af(2xi -a;2 -2:3) x 

{n^ + [Api + l)n + Ap\ + 2pi} - l{l + 1) 
+2Mxi(s2-l), (14) 

7„ = 2M{xi-X2){n^ +Apin + Apl- s^). (15) 

Here, s = 2 for the gravitational perturbations and s = 1 
for the electromagnetic perturbations. Note that in the 
present study, we do not consider the scalar perturbations 
because Leaver's method cannot be directly applied to 
the scalar case. Comparing the expanded eigenfunction 



(|ll(l with equation @, we can see that the eigenfunc- 
tion Hll|l satisfies the QNM boundary condition JHJ if the 
power series in equation (|ll(l converges for X2 <x<xi. 
This convergence condition is equivalent to the condi- 
tion written in terms of continued fractions [25> | , which is 
given by 

= /3o-^^--- (16) 

Therefore, we solve this algebraic equation to obtain the 
quasinormal frequency. 

III. NUMERICAL RESULTS 

In this study, we are concerned with asymptotic be- 
haviors of high-order QNMs, namely quasinormal fre- 
quencies with large imaginary parts, in almost the ex- 
tremal limit, in which two horizons are quite near. We 
therefore consider only the case of riKi < 10~^. For 
these small values of ri/ti, the continued fractions \i<6\ 
converge very quickly and Leaver's method works quite 
well even when quasinormal frequencies have quite large 
imaginary parts. Note that for moderate values of riKi, 
however, the convergence of the continued fractions gets 
worse and Leaver's method is applicable only for modes 
with smaller imaginary frequencies. For those cases, thus, 
some other techniques such as NoUert's method or a 
phase integral method [2^ should be used to obtain high- 
order QNMs. 

In order to check our numerical code, we have calcu- 
lated fundamental frequencies of the QNMs for several 
values of ri^i and have fitted the mode frequencies as a 
function of riKi with the polynomials defined by 

Re(a;ri) = riKi6o(l - &i?'i'«i) , 

Im(a;ri) = riKiCo(l - ciriKi) . (17) 

Recently, Cardoso and Lcmos and Maasse van den 
Brink [Tg analytically obtained the expansion coeffi- 
cients, which are given by 

f Jl{l + 1) - i for s = 0, 1 
[ + fors = 2, 



co = n+ 5i = ci - , (19) 

where n is the mode number. It is found that numer- 
ically obtained coefficients are in good agreement with 
coefficients given by (jl8|l and (|19|l . Note that our numer- 
ical results for the mode frequencies are consistent with 
those obtained by Moss and Norman [23, who studied 
low-order QNMs of the gravitational perturbations for 
full range of the parameter ri^i. 

First, let us discuss properties of the quasinormal fre- 
quencies for the low-order modes. In Figure 1, we show 



4 



the real parts of the frequencies for the low-order QNMs 
versus the imaginary parts of the frequencies for the 
gravitational perturbations with I ~ 2. In this figure, 
non-dimensional frequencies Re(cj/Ki) have been plot- 
ted as a function of Im(a;/Ki) for riKi = 10~^ and 
riKi = 5 X 10^'^, and the dashed curve indicates the 
approximate frequency for riKi — s- 0, derived by Cardoso 
and Lemos (see, also JJj). For low-order modes. 
Figure 1 illustrates how the frequencies of the QNMs be- 
haves when the mode number and/or the value of riKi 
is altered. Basic properties of the low-order QNMs are 
summarized as follows: For the modes associated with 
a small mode number, the real parts of the frequencies 
decrease with the increase of the imaginary parts of the 
frequency, even though the real parts of the frequency 
are constant in the approximation formula H17|l - H19() . In 
other words, the analytical approximation formula in the 
extremal limit (|17|I - H19(I is quite good for the fundamen- 
tal modes as long as riKi is small enough, while, as ex- 
pected in 10], this approximation formula gets worse as 
the mode number is increased even for small values of 
riKi. This means that formula (|17(l - (|19|l does not give a 
correct asymptotic value of the QNM frequencies in the 
limit of large imaginary frequencies. Similar properties 
can be seen for other perturbations having different s 
and I. In Figure 2, we show the same results as those in 
Figure 1 but for the electromagnetic perturbations hav- 
ing 1 = 1. It is observed that similar behaviors of the 
QNM frequencies are seen in the case of electromagnetic 
perturbations, too. 

Now, we explain our numerical results for the asymp- 
totic behavior of the QNM of nearly extremal SdS space- 
times in the limit of large imaginary frequencies. In Fig- 
ure 3, we plot the imaginary parts of the non-dimensional 
QNM frequencies, w/ki, of the gravitational perturba- 
tions with / = 2 as a function of the mode number n. 
In this figure, the model parameter is tiKi = 10^'^. The 
figure shows that an asymptotic form of Im(w/Ki) ^ n, 
which is similar to the analytical formula H17|) - (|19|l . is a 
good approximation for the imaginary parts in the limit 
of large imaginary frequencies. The same asymptotic 
form is inferred in all other QNM's we have calculated in 
the present study, regardless of the values of tiKi, /, and 
s. 

Let us next focus on the behaviors of the real part 
of the QNM frequencies in the limit of highly damped 
modes. In Figures 4 and 5, the real parts of the non- 
dimensional mode frequencies, lo/ki are plotted as a 
function of the imaginary parts of the frequencies up 
to sufhciently high-order modes for the Z = 2, 3 grav- 
itational and for the Z = 1, 2 electromagnetic pertur- 
bations, respectively. The results for the model param- 
eter of TiKi = 10"'^ are shown in both figures. It is 
found that the real parts of the frequencies show oscillat- 
ing behaviors as the imaginary parts of the frequencies 
are increased. It is important to note that similar os- 
cillating behaviors have been observed in the QNMs of 
a Reissner-Nordstrom black hole. (The quasinormal fre- 



quencies with large imaginary frequency of a Reissner- 
Nordstrom black hole can be given in terms of a solution 
of the algebraic equation 0, |3 01 1 

e'^" + 2 + 3e'^'^" = , (20) 

where (3 and k are contents determined with the mass 
and charge of the black hole (see, 0). As shown 
by Neitzke and Anderson and Howls H^, equation 
(|20|1 has an infinite number of solutions and some solu- 
tions of equation (|20|l show periodicity.) The behavior 
of the amplitude of the oscillating Re(a;/fvi) as a func- 
tion of Im(w/Ki) resembles that of QNM frequencies of 
a Schwarzschild black hole. For the gravitational per- 
turbations, the amplitude decreases for small values of 
Im(a;/Ki), approaches the imaginary axis of the complex 
frequency plan, increases again, and finally approaches 
some constant value. The asymptotic value of the am- 
plitude in the limit of highly damped modes seems to be 
non-zero constant, which is inferred as ~ 0.4. Therefore 
the limit of Re(w/Ki) as Im(w/Ki) — > oo appears not to 
exist for the gravitational perturbations. For the electro- 
magnetic perturbations, on the other hand, the ampli- 
tude of oscillating Re(tt'/Ki) decreases monotonically as 
Im(a;/Ki) is increased. Its asymptotic value in the limit 
of highly damped modes seems to be zero. This means 
that for the electromagnetic perturbations, the limit of 
Re(a'/Ki) as Im(w/Ki) goes to infinity seems to exist and 
to be zero. It is important to note that in a nearly ex- 
tremal SdS black hole, the asymptotic behaviors of the 
quasinormal frequencies in the limit of highly damped 
modes are independent of the angular quantum number 
I of the perturbations. Although we do not show the 
results for other values of riKi and the asymptotic be- 
havior of the QNM frequencies is not highly dependent 
on these parameters. In summary, our numerical results 
suggest that for the gravitational perturbations, the real 
parts of the QNM frequencies of nearly extremal SdS 
spacetimes do not go to any constant value in the limit 
of large imaginary frequencies because they show oscil- 
lating behaviors in the limit. For the electromagnetic 
perturbations, on the other hand, the real parts of the 
QNM frequencies seem to go to zero in the limit of large 
imaginary frequencies. 

IV. CONCLUSIONS 

We have calculated the high-order QNMs with large 
imaginary frequencies for the electromagnetic and gravi- 
tational perturbations in nearly extremal SdS spacetimes 
using Leaver's continued fraction method pol |. Our re- 
sults show that for low-order QNMs, analytical formulas 
in the extremal limit derived by Cardoso and Lemos [T^ 
and Maassen van den Brink |l7l | is a quite good approx- 
imation for the QNM frequencies as long as the model 
parameter tiKi is small enough. For high-order QNMs, 
whose imaginary frequencies are sufficiently large, on the 
other hand, this formula becomes inaccurate even for 
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small values of tiKi. Therefore, the approximation de- 
rived by Cardoso and Lemos cannot give correct asymp- 
totic behaviors of the QNMs in the limit of large imagi- 
nary frequencies (see, also T^). We also found that the 
real parts of the quasinormal frequencies have oscillating 
behaviors in the limit of highly damped modes. (Similar 
behaviors have been found in the quasinormal frequen- 
cies in a Reissner-Nordstrom black hole 10, 11, 12, 1^) 
The amplitude of oscillating Re(a;) approaches a non-zero 
constant value for the gravitational perturbations and 
zero for the electromagnetic perturbations in the limit of 
highly damped modes, regardless of values of I and riKi. 
This means that for the gravitational perturbations, the 
real parts of the quasinormal frequencies of nearly ex- 
tremal SdS spacetimes appear not to go to any constant 
value in the limit of highly damped modes. Therefore our 
numerical results suggest that Hod's conjecture is not ap- 
plicable for nearly extremal SdS black holes because the 
the limit of Re(Li;) as lm{Lu) —^ oo does not exist. 

Although we computed high-order QNMs whose damp- 
ing rates are quite large, all the QNMs we obtained in this 
study are still associated with a finite mode number but 
not infinity, because we investigated the properties of the 
QNMs with straightforward numerical approach. Thus, 
we cannot exclude the possibility that our numerical re- 
sults do not show correct asymptotic behaviors. Other 
approaches to examine asymptotic behaviors of QNMs 
in the highly damping limit are necessary, in order to 
confirm our results of the asymptotic behaviors. As for 
high-order QNMs with large imaginary frequencies for 
moderate values of ri^i. Leaver's method cannot be ap- 
plied straightforwardly. Therefore, other numerical tech- 
niques are needed to obtain QNMs with large imaginary 
frequencies. 
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APPENDIX: SUPERPARTNER AND 
ISO-SPECTRAL RELATIONSHIP BETWEEN 
AXIAL AND POLAR PERTURBATIONS IN SDS 
SPACETIME 

Non-radial gravitational perturbations of a SdS space- 
time obey a Schrodinger-type wave equation, given by 



dr'i 



where <^(+) and y(+) (0(-) and V^")) are the gauge in- 
variant perturbation function and the effective potential 
for polar (axial) parity perturbations, respectively. The 
effective potentials are given by 



2A 



r5(cr -I- 3M)2 

OA 

V'^-^ = [{c+l)r - 6M] , 



(A.2) 



(A.3) 



where c — {I + 2){l — l)/2 (for detailed derivations of 
the master equation IjA.ll) . see, e.g., |22,|3|)- shown 
first by Cardoso and Lemos 24], two potentials and 
are simply related through the relation, given by 



where (3 = 6M, n = Ac{c + 1), and 



/ 



2r3(cr + 3Af) 



(A.4) 



(A.5) 



This relation between two potentials are called the su- 
perpartner relationship. By virtue of the superpartner 
relationship, the perturbation function ('/''■"■*) can 

be written in term of 0'"' {4>'^^^) ^-nd its first derivative 
[23, given by 



(k ± 2iw/3)0(±) = (k- 



2/32/) 0(T)± 2/3^. (A.6) 



It is worthwhile to note that the superpartner relation- 
ship (|A.4p and l|A.5(l and the relations between two func- 
tions (jr^ (|A.6|I in a SdS spacetime have the same func- 
tional form as those in a Schwarzschild spacetime ex- 
cept for the definition of the function A Since 
A = A(r — ri)(r2 — r)(r — r3)r/3, the asymptotic form of 
the function / in the limit of — > ±00 is given by 



A(r2— ri)(ri— ra) 
6rJ(cri+3M) 

A(r2-ri)(r2-r3) 
6r^(cr2+3M) 



g2|Ki |r, 

„-2|K2k* 



as T ^ 



as 



(A.7) 



(±) 



0, 



(A.l) 



Then, it is easy to see that the function / has three prop- 
erties; i) smooth for —00 < < 00, ii) / and its deriva- 
tives of all orders vanish as ±00, iii) an integral 
1^00 Z*^^* exists. If the function / appearing in l|A.4|l 
satisfies three conditions above, as shown in p9l |. two po- 
tentials V^-^^ give the same transmission amplitude and 
the same quasinormal frequencies. In a SdS spacetime, 
therefore, axial and polar perturbations yield the same 
set of quasinormal mode frequencies. This iso-spectral 
properties in a SdS spacetime is attributed to the fact 
that gravitational perturbations associated with a spin 
s = — 2 in a SdS spacetime can be described with a sin- 
gle Weyl scalar *4 UM- 
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In a Schwarzschild-anti-de Sitter spacetime, exactly 
the same relations between polar and axial perturba- 
tions l|A.4|l - l)A.6|l obviously hold Yet, there is no 
iso-spectral property between polar and axial perturba- 
tions in a Schwarzschild-anti-de Sitter spacetime. In a 
Schwarzschild-anti-de Sitter spacetime, the master equa- 
tion IjA.ljl does not have an asymptotic solution given 
by e^*'^''* as r — > cx) and, furthermore, r» has a finite 
range. The boundary condition at spatial infinity there- 
fore must be modified. One of the plausible boundary 



conditions is that perturbation functions vanish at spa- 
tial infinity, even though there are other options for the 
boundary conditions [2^. If this boundary condition 
is taken at spatial infinity, in general, the transforma- 
tion (|A.6|I between polar and axial perturbations can- 
not hold this boundary condition. Therefore, the set of 
the quasinormal frequencies of polar perturbations in a 
Schwarzschild-anti-de Sitter spacetime is not the same as 
that of axial perturbations. 



[1] C. V. Vishveshwara, Nature 227, 936 (1970). 

[2] K. D. Kokkotas and B. G. Schmidt, Living Rev. Relativ- 
ity 2, 2 (1999). 

[3] H. P. Nollert, Class. Quant. Grav. 16, R159 (1999). 

[4] J. D. Bekenstein and V. F. Mukhanov, Phys. Lett B 360, 
7 (1995). 

[5] J. D. Bekenstein, in Proceedings of the Eight Marcel 
Grossmann Meeting on General Relativity, eds. T. Pi- 
ran and R. RufRni (World Scientific, Singapore 1999), 
pp. 92-111. 

[6] S. Hod, Phys. Rev. Lett. 81, 4293 (1998). 
[7] O. Dreyer, Phys. Rev. Lett. 90, 081301 (2003). 
[8] L. Motl, Adv. Theor. Math. Phys. 6, 1135 (2003). 
[9] H.-P. Nollert, Phys. Rev. D 47, 5253 (1993). 
[10] L. Motl and A. Neitzke, Adv. Theor. Math. Phys. 7, 307 
(2003). 

[11] E. Berti and K. D. Kokkotas, Phys. Rev. D 68, 044027 

(2003). ' 

[12] A. Neitzke, preprint | |hep-th/0304080t ^ 

[13] N. Andersson and C. J. Howls, preprint fgr-qc/03070201. 
[14] E. Berti, V. Cardoso, K. D. Kokkotas, and H. Onozawa, 

Phys. Rev. D 68, 124018 (2003). 
[15] S. Hod, Phys. Rev. D 67, 081501 (2003). 
[16] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 67, 084020 



(2003). 

[17] A. Maassen van den Brink, Phys. Rev. D 68, 047501 
(2003). 

[18] R. A. Konoplya, Phys. Rev. D 66, 044009 (2002). 

[19] V. Cardoso, R. A. Konoplya, and J. P. S. Lemos, Phys. 

Rev. D 68, 044024 (2003). 
[20] E. W. Leaver, Proc. R. Soc. London A402, 285 (1985). 
[21] I. G. Moss and J. P. Norman, Class. Quant. Grav. 19, 

2323 (2002). 

[22] F. Mellor and I. Moss, Phys. Rev. D 41, 403 (1990); 
[23] H. Otsuki and T. Futamase, Prog. Theor. Phys. 85, 771 
(1991); 

[24] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 64, 084017 
(2001). 

[25] W. Gautschi, SIAM Rev. 9, 24 (1967). 

[26] N. Andersson, Proc. R. Soc. London A439, 47 (1992). 

[27] V. Cardoso, J. P. S. Lemos, and S. Yoshida, Phys. Rev. 

D, in press (gr-qc/0309112 ); V. Cardoso, J. P. S. Lemos, 

and S. Yoshida, JHEP 0312, 041 (2003). 
[28] D. Birmingham, Phys. Lett. B 569, 199 (2003). 
[29] S. Chandrasekhar, The Mathematical Theory of Black 

holes (Oxford University, New York, 1983). 



7 



T-i-i-p-i-i-i-p-i-i-r 



r,K,=5xlO-3 



50 100 150 200 250 
Im(cj//cJ 



FIG. 1: Real parts of the non-dimensional QNM frequen- 
cies, u)/k\, given as a function of the imaginary parts of 
the frequencies for the 1=2 gravitational perturbations. The 
model parameters riKi are taken to be riKi = 10~* and 
riKi = 5 X 10~^. The frequencies obtained with the approx- 
imation formula in the limit of ri/ti — » are also shown as 
the dashed curve. 
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FIG. 2: Same as Figure 1 but for the I = 1 electromagnetic 
perturbations. 




FIG. 3: Imaginary parts of the non-dimensional QNM fre- 
quencies, uj/ni, given as a function of the mode number, n, 
for the gravitational perturbations associated with I = 2. The 
model parameter n/ti is taken to be riKi = 10~^. 
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FIG. 4: Real parts of the non-dimensional QNM frequencies, 
u)/k\, given as a function of the imaginary parts of the fre- 
quencies for the gravitational perturbations having / = 2 and 
I = 3. The frequencies of the QNMs associated with different 
I, I = 2, 3, are shown in each panel. The model parameters 
riKi is taken to be riKi = 10~^. 
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FIG. 5: Same as Figure 3 but for the electromagnetic pertur- 
bations having / = 1 and / = 2. 



